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Aims. Seeing-limited resolution in large telescopes working over wide wavelength range depends substantially on the turbulence outer 
scale and cannot be adequately described by one "seeing" value. We attempt to clarify frequent confusions on this matter. 
Methods. We study the effects of finite turbulence outer scale and partial adaptive corrections by means of analytical calculations and 
numerical simulations. 

Results. If a von Karman turbulence model is adopted, a simple approximate formula captures the dependence of atmospheric long- 
exposure resolution on the outer scale over the entire practically interesting range of telescope diameters and wavelengths. In the 
infrared (IR), the difference with the standard Kolmogorov seeing formula can exceed a factor of two. We find that low-order adaptive 
turbulence correction produces residual wave-fronts with effectively small outer scale, so even very low compensation order leads to 
a substantial improvement in resolution over seeing, compared to the standard theory. 

Conclusions. Seeing-limited resolution of large telescopes, especially in the IR, is currently under-estimated by not accounting for the 
outer scale. On the other hand, adaptive-optics systems designed for diffraction-limited imaging in the IR can improve the resolution 
in the visible by as much as two times. 
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1. Introduction 

Image blur of astronomical objects caused by terrestrial atmo- 
sphere is traditionally called "seeing". In the 2-nd half of the 
20-th century this ph enomenon has been understood and quan- 
tified (lYound 11974 . This understanding was based on con- 
sidering the distorted wave-fronts as a random stationary pro- 
cess wi th a power-law spectrum - the Kolmogorov-Obukhov 
model dTatars kii 1961; Roddier 1981). This theory describes the 
shape of the atmospheric long-exposure Point Spread Function 
(PSF) and many other phen omena by a single parameter, e.g. the 
Fried's coherence radius ro (Fried 1966). The theory predicts de- 
pendence of the PSF Full- Width at Half maximum (FWHM) eq 
on wavelength A and ro: 



The physics of turbulence implies that the spatial power 
spectral density (PSD) of phase distortions W^(f) (f is the spa- 
tial frequency in m ') deviates from the pure power law at 
low frequencies. A popular von Karman (vK) turbulence model 
( lTatarskiill96UIZiad et al.]|2000l:IConanll2000h introduces an ad- 
ditional parameter, the outer scale Lq: 



W^(f) = 0.0229 r^^'^ 



(If I' + Lfr 
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(2) 



£() = 0.976 A/ro. 



(1) 



In this paper we assume that ro and eq refer to observations at 
zenith. By adopting a standard wavelength A - 500 nm, we can 
replace ro with eo and this single parameter is nowadays usu- 
ally called "seeing". Here we use the term seeing in this precise 
sense, meaning eo at 500 nm at zenith. 

The success of this theory led most people to believe that 
the atmospheric parameters ro or eo actually exist and can be 
measured with high accuracy, given adequate means. In fact the 
match between real physical quantities like PSF or various sta- 
tistical estimates of distorted wave-fronts to the Kolmogorov- 
Obukhov theory varies from very good to poor, but it is never 
perfect. The concept of seeing becomes questionable if we push 
it too far. 
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Equation|2]is the definition of Lq. The Kolmogorov model cor- 
responds to Lo = oo. In the vK model, ro describes the high- 
frequency asymptotic behavior of the spectrum, and thus loses 
its sense of an equ ivalen t wavef ront coherence diameter as de- 
fined originally by iFriedl (Il966l) . Obviously, Eq. [T]is no longer 
valid as well. 

It remains an open question whether wave-front statistics ac- 
tually correspond to Eq.|2] Proving the vK model experimentally 
would be a difficult and eventually futile goal because large-scale 
wavefront perturbations are anything but stationary. However, 
it is firmly establishe d that the phase spectrum does deviat e 
from the power law ( IZiad etalj|2000l: iTokovinin et al]|2007l) . 
The Eq. |2] with additional parameter Lq provides a useful first- 
order description of this behavior Existing experimental data on 
Lo are interpreted here in this sense. 

In this paper, we study the modifications of Eq. [1] implied 
by the finite outer scale. Our analytical calculations are con- 
firmed by extensive numerical simulations. We show that for fi- 
nite L() the atmospheric FWHM becomes smaller than predicted 
by Eq. [T] and that this difference can be substantial. The practi- 
cal consequences for operation of large telescopes are discussed. 
The lack of low-frequency power is typical not only for the vK 
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turbulence, but also for partially corrected wave-fronts result- 
ing, e.g., from tip-tilt compensation (fast guiding) or low-order 
adaptive-optics (AO) correction. Such a correction leads to small 
effective Lq. We apply the same analytical treatment to this case 
and study the shrinking of the PSF halo under partial AO com- 
pensation. 

2. Analytical treatment 

The calculation of the long-exposure PSF is done by multiplying 
the telescope optical transfer function (OTP) by an additional 
term, the atmospheric OTF: 



T„(u) = exp[-0.5D^(^u)], 



(3) 



where u is the angular spatial frequency (in inverse radians), A is 
the imaging wavelength, and Z)^(r) is the phase structure func- 
tion (SF) (Goodman 1985; Roddier 1981). This expression is 
general, applicable to any turbulence spectrum and any telescope 
diameter. In the case of a large ideal telescope with diameter 
D » ro the diffraction can be neglected and the long-exposure 
OTF and PSF are accurately described by Eq. [3] 

The analytic expression for the phase stru cture function 
in the von Karman model can be found in (IConanI \m)dt. 
IConsortini et al.|[T972l : lTokovininl2002l) . For infinite Lq, it trans- 
forms into D^(r) = 6.88(r/ro)^/l 

Figure[T](top) plots the SFs for Kolmogorov and vK with the 
same r^. In the latter case the SF saturates at r > Lq, reaching 
asymptotically the level 0. 17(Lo/'"o)^^^- It reaches half-saturation 
at r = 0.17Lo- The Kolmogorov SF with the same ro crosses the 
vK saturation level at r = 0.109Lo. This tells us that the effect 
of finite outer scale is strong at distances much shorter than Lq, 
and that it would be misleading to compare Lq directly with the 
telescope diameter. 

Putting the vK SF into Eq.[3l we find that for finite Lq Ta(u) 
does not go to zero at large arguments, therefore its inverse 
Fourier transform (the PSF) formally does not exist. However, in 
the case when ro <s Lq this level is small and it can be neglected. 
In Fig. [T] (bottom) we compare the PSF profiles for different val- 
ues of L()/r(), including Lq - oo. 

A first-order approximation of the FWHM of atmospheric 
PSFs (fivK) und er vK turbulence has been suggested by 
iTokovininI (120021) : 



evK ^ so Vl -2.183 (ro/Lo)' 



,0.356 



(4) 



This formula is valid for Lo/ro > 20 to an accuracy of +1%. We 
remind the reader that while ro depends on the wavelength, Lq 
does not. At smaller Lo/ro values, the atmospheric PSF develops 
a strong core-halo structure, and its FWHM becomes less and 
less meaningful. The actual PSF in a telescope is a convolution 
of the atmospheric blur with diffraction, aberrations, guiding er- 
rors, etc. As neither of these factors is described by a Gaussian, 
calculation of the combined FWHM as a quadratic sum of indi- 
vidual contributions is not accurate. 

Similarly a formula for the FWHE, ha lf-energy diamet er (J3yk), 
can be derived with the same accuracy (lTokovimnl2002l) : 



ySvK ^ySo Vl- 1-534 (roW^^^- 



(5) 



where jSq is the diameter of the circle that contains one-half of 
the total PSF energy in the Kolmogorov theory (J3q = 1.15/l/ro). 

The following section gathers results obtained with exten- 
sive numerical simulations to confirm the reliability and vaUdity 
domain of Eq.Hl and thereby Eq.|5] 
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Fig. 1. Top panel: comparison of the von Karman (full line) and 
Kolmogorov (dashed line) phase structure functions with the 
same ro. Bottom panel: normalized atmospheric PSFs for dif- 
ferent Lo/ro ratios. 



3. Numerical simulations 

3. 1 . Random wavefronts 

The atmospheric turbulence is generated with 1000 uncorrelated 
phase screens on 8192 x 8192 array equivalent to a 100 me- 
ters width physical size (pixel size 12.2 mm). The principle of 
the generation of a phase screen is based on the Fourier ap- 
proach (McGlamery 1976) : randomized white noise maps are 
colored in the Fourier space by the turbulence PSD (Eq.|2]i; the 
inverse Fourier transform of an outcome corresponds to a phase 
screen realization. The large size of the simulated phase screens 
is mandatory to correctly sample the Lo and to compute PSF for 
large telescopes. The simulations consider several Lq cases (10, 
22, 32.5, 50, 65 m and oo). 

Several investigations have been carried out on the phase 
screens to ascertain that their statistics correspond indeed to the 
input parameters ro and Lo- For example, we compared the phase 
variance, and the variances of the first 100 Zernike coefficients 
for D = 42 m ( ro=12.12cm a nd Lo - 65 m) with their expected 
values given bv lConad (120001) and found a good agreement. The 
phase variance matches expectation to within 1.7%, while the 
variance of tip and tilt components meets their theoretical values 
to within 1.3 and 0.7% accuracy, respectively. The case Lo = oo 
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Fig. 2. The atmospheric FWHM of simulated long-exposure PSFs 0.^ versus telescope diameter for several Lq values (A=Q.5fim, 
SQ = 0.83"). The blue curves trace the diffraction FWHM 0dif = ^/D. 



is particular: the variance of the tip and tilt coefficients does not 
fit their theoretical values, corresponding instead to a finite outer 
scale in the range of 200 m to 500 m. This is a consequence of 
the finite size of the simulated phase screens. 

Several telescope diameters have been considered ranging 
from 10 cm to 42 m. The wavelength domain ranges from the 
U-band to M-band, while the seeing ranges from 0.1" to 1.8". 

All simulations involve Fast Fourier Transforms (FFT) of 
8192 X 8192 aiTays to generate the long-exposure PSFs (over 
1000 realizations). The same sets of phase screens is used for 
all telescope diameters. As a result of the very large arrays in- 
volved in order to handle both phase screen statistics and alias- 
ing effect (e.g. large-aperture cases), small telescope diameters 
(< 1 m) may suffer from coarse pupil sampling. The effect of 
speckle structure is also stronger for small D, causing a larger 
random scatter in the results. 

3.2. Measurement of the FWHM 

We determined the FWHM of the simulated long-exposure PSFs 
Opsp in the following way. The PSFs were first azimuthally aver- 
aged. The 10-th order polynomial was then fitted to this curve, 
and the radius where it crosses the 1/2 of the maximum intensity 
was determined. The outc omes of th is routine has been com- 
pared to another algorithm (lKolbll2005i) applied on the same set 



of PSFs, and both gave similar values (e.g. +1 pixel for D = 8m 
and Lo = 22m, i.e. 0.006"). 

The simulated PSFs are broadened by diffraction and thus 
are not directly comparable to Eq. |4] We approximately account 
for this by subtracting quadratically the diffraction FWHM fljif = 
A/D, 

0. « V^psf' - 0dif'. (6) 

This gives a good approximation to EvK as long as the diffrac- 
tion blur is small, D s> ro, but fails at small diameters, as men- 
tioned above, because the individual broadening factors are not 
Gaussian. This explains why our results for small D are inaccu- 
rate. 

3.3. Outer scale and telescope diameter 

The first series of simulations aims at defining the general trend 
of atmospheric FWHM in large telescopes in the presence of 
finite outer scale. We compare 9-^ to Eq. |4]and to the seeing sq, 
fixed at 0.83" in this case. Some results are presented in Fig.|2] 
From Fig. |2]it is straightforward to see that 0a < eo in all 
cases, even for Lq - oo because all simulated wave-fronts have 
finite outer scale. As expected, the validity of Eq.|4]is confirmed, 
except for the small D where our treatment of diffraction is too 
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Fig. 3. Dependence of 6a on wavelength (top, sq = 0.83") and 
seeing (bottom, A = O.Syum ). Other parameters: Lq - 22 m, 
D = 8m. 



crude. All these cases correspond Lo/ro > 80 where the effect of 
the finite Lq is still mild. 



3.4. Wavelength and seeing dependence 

For the second series of simulation, we have considered an 8- 
m telescope, a fixed outer scale Lq - 22 m, and 0.83" seeing 
at O.Sfim, while the imaging wavelength is varying from the U- 
band to the M-band (from 0.365 to 4.67 //m). The results are 
presented in Fig. [3] top. Note the stronger dependence of 6^ on 
wavelength, compared to the Kolmogorov case. The third series 
of simulation considers the same Lq and D, while the seeing con- 
ditions are evolving (Fig. [3] bottom). The agreement with Eq.|4] 
is demonstrated for both wavelength {Lo/ro > 10) and seeing 
dependence (Lo/ro > 20). 



3.5. Discussion 

The previous subsections gave general results for the atmo- 
spheric FWHM in the presence of a finite outer scale. In order 
to relate these results to actual situation, we discuss the partic- 
ular case of the 8-m Very Large Telescope at Paranal (Chile), 
and assuming standard seeing conditions (0.83" at 0.5 i-im), and 
median outer scale value (Lo - 22 m, i.e. Lo/ro - 180). Results 



Fig. 4. Top: Ratio between seeing so and FWHM (e, /j-) as a func- 
tion of the wavelength for several Lq cases. Bottom: similarly, ra- 
tio between Po and FWHE (fi^K). For both plots, seeing is 0.83" 
at 0.5yL(m. 



shown in Fig.[3]indicate that the FWHM of the von Karman PSF 
is reduced by 19% compared to standard theory {eq) in the vis- 
ible. It is even more dramatic in the near-IR, where the FWHM 
(syK) is reduced by 29.7% (H-band) and 36.3% (K-band). 

In the same way, Fig.|4](top) quantifies ratio of sq by s^k (i.e. 
Eq.[T]to Eq.|4]i but for several Lq values ranging from the Paranal 
median value (including the l-cr outer scale values, 1 3 and 37 m) 
to 150 m. The difference with the standard Kolmogorov seeing 
formula is substantial and can exceed a factor of two in the IR. 
Likewise, Fig. |4] (bottom) compares ratio of Po by /3vk- 

4. Resolution under partial compensation 

In analogy with the finite outer scale impact, we discuss here the 
consequences of a reduction of the low-frequency content of the 
phase perturbation spectrum generated by AO partial correction, 
or tip-tilt compensation. 

The purpose of Adaptive Optics (AO) systems is to compen- 
sate atmospheric wave-front distortions and to reach diffraction- 
limited resolution. To do this, the actuator spacing d (or an equiv- 
alent measure of AO compensation order) must be of the order 
of 2ro or smaller (Roddier 1998). However, AO systems which 
do not fulfill this condition still improve the resolution. A good 
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example is the use of the low-or der AO system PU EO for ob- 
servations at visible wavelengths (iRigaut et alJll998]) . A resolu- 
tion gain up to two times has been reported. To our knowledge, 
shrinking of the atmospheric PSF under partial compensation 
has not been explored in a systematic way. 

Residual wave-fronts after AO compensation contain high- 
frequency ripple, whereas the perturbations at spatial frequen- 
cies less than fc - 1 / {2d) are coiTected. This can be modeled 
by high-pass filtering of the atmospheric PSD. The form of this 
filter varies, depending on the AO system. The calculations here 
only illustrate the principle and should be repeated for each AO 
system if an exact result is sought. We model the AO compensa- 
tion by a multiplicative factor F: 

F{x)^xl{\+x), x^{\i\l0" (7) 

with m - 6. The PSD (Eq.|2|i is multiplied by F, the SF is cal- 
culated and used to compute the residual PSF in the same way 
as for the vK spectrum. The SF saturates at r » t/, reaching 
the value 2cr^ = 0.62(ii/ro)^^^. The shape of the SF and the sat- 
uration value depend on the filter F{x). We experimented with 
several filters and have chosen Eq. Q with m - 6 because it 
matches approximately the known formula cr^ - 0.35(t//ro)^^"', 
lRoddierl(ll998h . Comparing this to the saturation level of the vK 
SF, 0.17(Lo/'"o)^''', we may state that the effective outer scale of 
the residual wavefront is ~ 2d. 

Once the SF saturates at 2cr^, the atmospheric OTF reaches 

a constant level T^i^ = exp(-cr^). We can represent such an 
OTF as a sum of the constant term and a decreasing part. This 
corresponds to the sum of a diffraction-limited PSF scaled by 
S - exp(-cr^) and a wide residual halo. The shape of the halo 
can therefore be calculated by replacing the atmospheric OTF 
with (Ta - rmin)/(l - Tmin) out to the distance where the mini- 
mum Tmin is reached, and setting it to zero for larger frequencies. 
Note that we re-normalize the halo OTF to one at the coordinate 
origin. 

The FWHM Oao and the diameter of a circle containing half 
the energy (FWHE) have been computed for the halo of par- 
tially compensated PSFs. We compare these parameters to the 
non-compensated (Kolmogorov) PSFs in Fig.|5] Even when the 
actuator spacing d is much larger than ro and the AO system 
does not perform well in the classical sense (S ~ 0), the gain 
in FWHM and FWHE is already substantial. Maximum resolu- 
tion gain so/0Ao is reached at d/ro ~ 4, when the coherent PSF 
core is still very weak. As the compensation order increases fur- 
ther, increasing fraction of energy goes into the core, the PSF 
halo becomes weaker and wider At small d and high S , the halo 
becomes even wider than the un-compensated atmospheric PSF, 
being produced by residual phase errors at spatial scales smaller 
than ro. 

Tip-tilt correction is a particular case of low-order AO com- 
pensation. It is well k nown that maximum resolution gain is 
achieved at D/ro ~ 3.6 (lFriedlll96 6). The gain studied here does 
not depend on the telescope diameter D, but rather on the di- 
mensionless parameter d/ro, in full analogy with the effect of 
the outer scale. 

All three effects - outer scale, partial AO correction and 
tip-tilt compensation - reduce the low-frequency content of the 
phase perturbations spectrum. When they act together, the gain 
in resolution over Kolmogorov turbulence is not cumulative. For 
example, with finite outer scale the tip-tilt fluctuations become 
smaller and their correction achieves a smaller resolution gain. 
Similarly, the resolution gain from partial AO correction (Fig.|5]l 
in fact will be less because of the finite Lq- 
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Fig. 5. Gain in the FWHM (fuH line) and FWHE (dotted line) 
diameter of the PSF halo (compared to the Kolmogorov PSF) re- 
sulting from partial AO compensation (outer scale not included). 
The dotted line shows the coherent energy S . 



5. Conclusions and discussion 

This study is largely motivated by the confusion between see- 
ing and the FWHM of long-exposure images in large telescopes, 
also often called delivered image quality (DIQ). In an ideal large 
telescope (no aberrations, internal turbulence and wind shake), 
the DIQ is always less than predicted by the standard theory, 
owing to the finite turbulence outer scale. 

The seeing is usually m easure d by the Dif ferential Image 
Motion Monitors (Martin 119871; ISar azin & R oddier 1993 
DIMMs). This method is sensitive to small-scale wave-front dis- 
tortions an d provid es estimates of ro which are almost indepen- 
dent of Lo (IZiad et al. 1994)0. Using the standard theory, we will 
over-estimate the FWHM expected in a large telescope. As the 
PSF is broadened by non-atmospheric factors, this mismatch can 
hide telescope defects. This is particularly dangerous in the IR, 
where the difference with the standard theory is large. Therefore, 
reaching a truly seeing-limited telescope performance in the IR 
requires mandatory account of finite Lq. Stated in other words, 
our telescopes could perform better than we think they should 
based on the standard theory and DIMM measurements. 

If, on the other hand, we want to deduce atmospheric see- 
ing from the width of the long-exposure PSF, the situation is 
reversed. The actual seeing is worse than we think it is. The 
effect of finite Lq is apparent for all telescope diameters. A si- 
multaneous measurement of Lq is thus required to be accurate. 
Estimating seeing from the width of the spots in active-optics 
Shack-Hartmann sensor (long exposures) should be done with 
these circumstances in mind. 

As internal telescope defects and outer scale act in opposite 
directions, they can partially compensate each other An agree- 
ment between DIMM m easureme nts and DIQ can thus be found 
where it should not be (ISarazin & Roddier 1990). By compar- 
ing simultaneous PSFs at visible and in the mid-IR, it is possi- 
ble to extract two parameters, eq and Lq, assuming that the tele- 
scope's contrib ution to the image degradation can be neglected 
(iTokovinin et al.ll2007l) . 



' Seeing monitors based on the absolute image motion are affected 
by finite Lq, giving biased, larger ro values. 
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